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We explore the influence of an arbitrary external potential perturbation V on the spectral prop- 
erties of a weakly disordered conductor. In the framework of a statistical fleld theory of a nonlinear 
cr-model type we flnd, depending on the range and the proflle of the external perturbation, two 
qualitatively different universal regimes of parametric spectral statistics (i.e. cross-correlations be- 
tween the spectra of Hamiltonians H and H -\-V). We identify the translational invariance of the 
correlations in the space of Hamiltonians as the key indicator of universality, and find the connection 
between the coordinate system in this space which makes the translational invariance manifest, and 
the physically measurable properties of the system. In particular, in the case of localized pertur- 
bations, the latter turn out to be the eigenphases of the scattering matrix for scattering off the 
perturbing potential V . They also have a purely statistical interpretation in terms of the moments 
of the level velocity distribution. Finally, on the basis of this analysis, a set of results obtained 
recently by the authors using random matrix theory methods is shown to be applicable to a much 
wider class of disordered and chaotic structures. 

PACS numbers: 05.45.Mt,73.21.-b,03.65.Sq 
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I. INTRODUCTION 

The statistical approach to the study of the quantum 
spectra of complex Hamiltonians H was pioneered in the 
1950's in connection with the study of resonances in nu- 
clear scattering Hill. The approach was formalized by 
Dyson who identified three universal types of spec- 
tral statistics as determined by the symmetry properties 
of the corresponding ensembles of random Hamiltonians. 
The theory of Wigner and Dyson, together with its ex- 
tensions, is usually referred to as random matrix the- 
ory (RMT). Within RMT, spectral correlations of each 
type (symmetry class) of statistics are described by the 
corresponding set of universal functions, once the ener- 
gies are rescaled by the average level spacing A. The 
latter is the only non-universal parameter of the the- 
ory, where non-universality implies dependence on the 
specifics (other than the symmetry class) of the distribu- 
tion functions defining the Hamiltonian ensemble. The 
symmetry classes are labeled by the index (3 which takes 
the values (3 = 1 for real (time-invariant) Hamiltonians 
(so-called orthogonal class), /3 = 2 for complex Hamilto- 
nians (unitary class), and /3 = 4 for Hamiltonians with 
broken spin-rotation symmetry (symplectic class). 

For example, in the unitary class (/3 = 2), the nor- 
malized two-point correlation function of the density of 
states (DoS), v{e) = tr^(e — H), 



where 



(1) 



can be expressed as 



i?2(s) = l + -Re / dXi dXe''''^^^'-^'^ 

= 5{s) + l~k^{s) 



sin(7rs) 



(3) 



(2) 



s = ri/A, and where 3+ = s + iO. The delta function 
in ^ represents the autocorrelation of an energy level 
with itself; the constant term is the disconnected con- 
tribution, and the oscillating term corresponds to non- 
trivial correlations between different levels. Similar, al- 
though slightly more convoluted, expressions are well 
known for other symmetry classes. 

As became clear in the subsequent development of the 
theory, Wigner-Dyson statistics possesses a remarkable 
universality: beyond the manifestly model constructions 
of RMT, Wigner-Dyson statistics occurs both in quan- 
tum weakly disordered systems 0, 0] = + ^^dis, 
where the statistical description follows naturally from 
the ensemble of realizations of the disordered part of 
the Hamiltonian i?dis, as well as in individual classi- 
cally chaotic Hamiltonians ^ upon averaging over a suf- 
ficiently wide spectral window. The status of the last two 
statements is somewhat different: while the former has 
been rigorously established '3, ^ , the latter statement is 
known as the Bohigas-Giannoni-Schmit (BGS) 5] con- 
jecture. Crucially, the existence of a statistical ensemble 
allows average properties of the disordered Hamiltonian 
to be formulated in the framework of a quantum field 
theory of nonlinear sigma model type (NLcrM) valid in 
the limit g ^ I. Here, the dimensionless conductance g 
denotes a non-universal parameter which depends on the 
deterministic part of the Hamiltonian Hq - and through 
it on the sample geometry - as well as on details of the 
distribution of i?dis- By contrast, properties of an indi- 
vidual system are much less amenable to treatment by 
means of a statistical field theory and, despite recent en- 
couraging progress i a convincing formal proof of the 
BGS conjecture is lacking. The existence of a proof of 
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Wigner-Dyson statistics in quantum disordered systems 
generates a dichotomy of approaches: if universahty is 
assumed or estabhshed by other means, RMT calcula- 
tions can be employed to obtain specific answers. On the 
other hand, NLcrM calculations often represent the only 
available route to establish universality and to explore 
deviations from it (typically as expansions in 1/g). 

A corollary of the above universality is the fact that, 
apart form a possible modification of the average DoS 
A~^, the spectral statistics of H = Hq + H^is and 
Hy — Ho + Hdis + V are indistinguishable unless 
V is so large that Hq + i?dis generate only perturba- 
tive corrections to its eigenvalues (the same condition, 
of course, applies to the relation between Hq and i?dis)- 
Such universality dictates that the statistics of spectral 
lines tracing the evolution of eigenvalues as functions of 
a coordinate X along a typical curve in the Hamiltonian 
space connecting Hq to Hq + V must possess some kind 
of stationarity j^j. This stationarity is reflected, in turn, 
in the universality of parametric spectral correlations 
(i.e. cross-correlations between the spectra of H and Hy 
averaged over the realizations of i?dis)- Indeed, as first 
pointed out in Refs. [Hill 113, 

for a broad class of V^s 
the parametric spectral correlators are expressed via an- 
other set of universal functions (again determined only 
by the symmetry of the Hamiltonian ensemble) which 
depend, in addition to the level density, on a single ex- 
tra non-universal parameter - the dispersion C (0) of the 
distribution of level velocities A~^dei/dX = A~^dx£i, 
where is a typical eigenlevel of H. (The use of the 
notation C(0) for the dispersion of the distribution of 
level velocities is explained by the fact that it is the 
limiting value of the level velocity correlation function 
C{X) = A-^{^xe^{X)^x£^{X + X)) If the 

curve is a 'straight line' H + XV, the rescaled parameter 

xo^Xy/CiO) (4) 

is an effective measure of the strength of XV. Eq. Q 
is, in essence, a relation between two phenomenological 
parameters, xq and C(0). It is supplemented by a 'mi- 
croscopic' definition of xq in terms of XV 0, 0| 

xl^^X'iiVUV)) , (5) 

where the double angle brackets denote the matrix ele- 
ment of n on the vector space spanned by V, and the 
definition of the (non-universal) operator 11 strongly de- 
pends on Hq and the distribution of -ffdis- In a disor- 
dered metal, 11 becomes a resolvent of the diffusion oper- 
ator 0,^3, ^'^'i' if ^ is diagonal in the coordinate rep- 
resentation, the matrix element simplifies to {{VHV)) — 
L-^'^ J drdr'y(r)n(r,r')F(r'), where n(r,r') is the ma- 
trix element of 11 in the coordinate basis, L is the linear 
size of the sample, and d are the space dimensions. The 
parameter xq can be thought of as a scalar norm (in the 
space of Hamiltonians) of the 'distance' between H and 



H + XV, and Eq. lO can correspondingly be viewed as 
fixing a particular prescription for the norm. 

The importance of Eq. is underscored by the fact 
that in practice, measuring V directly is often impossible 
or very difficult, and the Green function of the diffusion 
operator is itself a complicated object in an irregularly 
shaped sample. Consequently, direct evaluation of jSJ is 
seldom feasible. On the other hand, collecting statistics 
on level velocities is relatively straightforward, so that 
Eq. Q is the only link between the formal parameter xq 
and an experimentally measurable quantity. 

[A brief remark on notation: we denote the paramet- 
ric perturbation as XV whenever the emphasis is on the 
continuous evolution of the spectrum of H -I- XV as a 
function of X along a 'straight line' < X < 1. In 
what follows non-linear parameter dependence of V will 
become important, so it will be convenient to revert to 
denoting the perturbation as V with its parameter de- 
pendence implicit from the context.] 

The 'broad class of V's' which was referred to above 
is distinguished by the following properties. First of all, 
the main feature of the class of V's to which the theory 
of Refs. [lol [l3 applies is that the level velocity distri- 
bution is Gaussian, and is thus fully characterized by its 
second cumulant. It is, of course, generic in some sense, 
as the Gaussianity of the level velocity distribution is the 
result of the central limit theorem which comes into force 
for any V which is global, or extended, i.e. when V is a, 
generic 'full' matrix in the Hilbert space defined by a 
typical realization of Hq + H^is (a more formal definition 
is tr V^''^/(tr F^)^ in the thermodynamic limit). In 
application to disordered metals, such perturbations are 
easily realized as, e. g., inhomogeneous electric or mag- 
netic fields acting on the whole, or a substantial part of, 
the sample volume. 

Secondly, the fact that xq is linear in X implies that V 
is, in some sense, small. More precisely, individual ma- 
trix elements of V are small as compared to the mean 
level spacing A, even as the overall effect of V is finite 
due to its global extent. It cannot be too strongly em- 
phasized that, for a global perturbation, the requirement 
that its matrix elements are small is not an artificial con- 
straint, but a self-consistent condition for the existence of 
finite correlations between the spectra of H and H + V: 
if this condition is violated, the spectrum of H + V is 
so 'scrambled' with respect to the spectrum of H that 
the residual correlations vanish as inverse powers of the 
system volume. 

There exists, however, a class of perturbations V which 
don't quite fall into the paradigm of Refs. 0, ^^]. In 
many practical applications, for example, V cannot be 
characterized as extended. A bistable defect jumping 
between two nearby configurations, a local STM tip, de- 
fects created under irradiation, etc., belong to a very dif- 
ferent type of a perturbing Hamiltonian. On an intu- 
itive level they are easily seen to be 'local', and formally 
tr T^^/(tr y^)^ is finite for such perturbations. Even more 
crucially, the corresponding distribution of level velocities 
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is not Gaussian [e.g., it is Poissonian for a moving local 
defect in a disordered metal for = 2, see Eq. be- 
low], so the results of Refs. 0,0] need to be modified. 

Parametric spectral correlations induced by non-exten- 
ded perturbing potentials have been studied previously 
in the RMT framework 0, 0, ^| , leading to a com- 
prehensive set of results for a large class of parametric 
correlation functions in the unitary ensemble. As noted 
above, the use of RMT presupposes universality. At the 
same time, matrix elements of a local perturbation need 
not be small. As a typical example, let us consider the 
potential arising in the x-ray emission/absorption spectra 
in disordered metals, and the associated issue of Ander- 
son orthogonality catastrophe . The potential U (r) of 
a charged ion created after an electron has been knocked 
out of an inner shell by a photon may typically have 
a magnitude comparable to the Fermi energy ep and 
a spatial extent of the order of the Fermi wavelength 
Af = 27r/pF. For comparison, an exarnple of a global po- 
tential to which analysis of Refs. 0, 0| fully applies is 
external magnetic field generating a unit of flux through 
the sample area. The corresponding effective potential 
is easily estimated as (e/mc)pA ~ ep/ppL, where A is 
the vector potential. In the former case, the potential, al- 
though spatially localized, is larger by a factor of square 
root of the sample volume (in d = 2 dimensions). The 
qualitative difference between the two types of perturba- 
tions is vividly illustrated in Figs. ^ and [5] which depict 
the evolution of a set of levels of a typical realization of 
Hq + i?dis + XV as a function of X for the cases of global 
and local perturbations, respectively. Note, for example, 
that in the case of a local perturbation, even very large 
values of X do not lead to a complete 'scrambling' of the 
spectrum. 

Given these qualitative differences, the universality 
of the corresponding parametric correlations cannot be 
automatically deduced from the analysis of Refs. 0, 
IT^ . In this paper, we utilize the NLctM approach to 
prove the universality of the results in 0, Q, ^| and 
to establish their limits of applicability and, crucially, 
their parametrization in terms of non-universal (i?o- 
dependent) quantities. We also establish a precise crite- 
rion for the crossover between the localized 0, and 
extended 0,0] regimes. The NLctM also affords an ex- 
tension of some of the previously obtained results to the 
orthogonal ensemble, although a complete counterpart of 
the results of Refs. [T^, 111, [l3 in the orthogonal case is 
unattainable at present. 
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FIG. 1: Parametric dependence of the energy levels of a tight- 
binding Anderson Hamiltonian on a 20 x 20 lattice as a func- 
tion of the strength X of a fixed global perturbation (arbitrary 
units). To place the system in the unitary symmetry class, we 
have imposed an external random vector potential in addition 
to the random scalar potential. 
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FIG. 2: Parametric dependence of the energy levels of a time- 
reversal-invariant tight-binding Anderson Hamiltonian on a 
20 X 20 lattice as a function of the overall strength X of an 
external perturbation localized on two fixed sites (arbitrary 
units). 



II. PARAMETRIC SPECTRAL CORRELATION 
FUNCTIONS 

The parametric analog of Eq. - and the basic object 
to which we devote most of the attention in this paper - 
is the cross-correlation function 

Rn{n, V) = h^{v{e + 17/2; 0) v{e ~ n/2; V)) , (6) 



where i^(e; V) = tr5(e — H — V). The notation is dic- 
tated by the fact that Ru is the special case of the 
multi-point correlation function Rnm involving n densi- 
ties ^(£i;0) and m 'shifted' densities v{£j]V). (Within 
this classification scheme, i?2 is more properly denoted 
as -R20-) Its universal form (in the case of random Hamil- 
tonians of unitary symmetry) is given, as established in 
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Refs. [13,112 for global V, by the following expression: 
Riiis,V)-l 

-, poo pi 

= -Rey d\i J dAe*'''+(^^-^)-''(^'^^'^'') , (7) 

where 

ct(A, Ai; a;o) = (7gi(A, Ai; xo) = TT^xUXj - A2)/2 , (8) 

and xq is given by Eq. Q . Note that the resolvent of the 
diffusion operator n(r,r') is a long-range object, so that 
the extended character of V is essential in the structure 
of Eq. ||SJ|. Note also that 11 is smooth on the scale of 
the mean free path £, suppressing the contribution of the 
fast components of to lO- 

The basic conclusion of the analysis presented in this 
paper is that in the case of a local perturbation (or a per- 
turbation which has both global and local components) 
the overall structure of Eq. (TJl is preserved, while the 
bmction a(X, Xi; xq) is generalized to 

(t(A, Ai;a;o,x) = a-gi(A, Ai; xq) + crioc(A, Ai; x) 

r 

CTioc(A, Ai; x) = ^ [In (1 + iXaXi) - In (1 + ix^X)] , 

a=l 

(9) 

where the 'local' part aioc of cr depends on a vector of 
parameters x. The vector x has a finite number r of 
components, reflecting the finite extent of the local com- 
ponent of V. We also find an additional contribution to 
Xq, which properly accounts for the fast spatial oscilla- 
tions in the extended parts of V, 

6x1 = 4i^^ J drdr'V{r)fj{r - r')Vir')X^ , (10) 

where D — is the average local DoS in a d- 

dimensional sample. The Friedel function /d(r) repre- 
sents the averaged Green function of a disordered system, 
and is given by 

/d(r) = J<i/2-l(PFr)e-'■/2^ (1^) 

where J,, is the Bessel function of 77-th order. 

In what sense is Eq. © universal, depending, as it 
does, on a whole list of parameters? To answer this ques- 
tion one needs to clarify whether analogs of Eq. Q can 
be defined, i.e., whether these correlation functions are 
determined by an (at most finite) set of phenomenologi- 
cal parameters rather than the thermodynamically large 
number of variables in the pair of matrices -ffo and V. 

As will be explained shortly, Xa are nonlinear func- 
tions of XV. In other words, the transformation analo- 
gous to Q, if it exists at all, is non-linear - in contrast 
to the simple rescaling required in the global case. The 
nonlinearity means that, unlike the case of a global per- 
turbation, the overall magnitude of the perturbation X is 



no longer a 'natural' variable with respect to which level 
velocity should be defined. It is not even a priori clear 
that such a 'natural' variable exists. 

There is, however, an extra requirement that can be 
imposed to make the parametrization unique. To help 
motivate it, let us note that, in the global case, a as 
a function of xq is essentially a Taylor expansion where 
terms of the order higher than 2 vanish in the thermody- 
namic limit. As a result, correlation functions between 
H + XiV and H + X2V for arbitrary Xi and X2 are 
functions of the (properly rescaled) difference X2 — Xi. 
This translationally invariant parametrization is a direct 
consequence - in fact, the primary manifestation - of 
the principle of stationarity discussed in the Introduc- 
tion. On general grounds, such stationarity should char- 
acterize parametric correlations induced by local pertur- 
bations as well. Thus, there must exist a 'natural' set of 
variables y(x) possessing the following property: there is 
at least one curve connecting Hq and Hq + V such that 
the correlations between any two Hamiltonians on this 
curve corresponding to two different values Xi and X2 
are functions of the differences y(x2) — y(xi). The distri- 
bution of level velocities along this curve is independent 
of the position on the curve. In fact, for r > 1, such 
curves are not unique, but rather span an r-dimensional 
manifold. 

It is worth emphasizing that the existence of such 
parametrization is a non-trivial property even for a single 
parameter (r — 1). For example, let us consider an arbi- 
trary function f{Xi, X2). The existence of a 'natural' set 
of variables is equivalent to the requirement that for any 
Xi and X2 the function f{Xi,X2) can be represented 
as g{h{X2) — h{Xi)) for some functions g and h. Such 
a representation clearly exists only for a very restricted 
subset of all possible /. As follows from the analysis be- 
low, the 'natural' parametrization does indeed exist, and 
it is effected by the substitution 

Ua = arctanxo . (12) 

As will be shown in Section UTTl the eigenvalues of 

the reactance matrix 

n = i{S -l)iS + (13) 

where the scattering matrix S describes scattering off V 
with free propagation controlled by the Green function 
oi H = Hq + Hdis averaged over i?dis- The emergence 
of the reactance matrix as a key object governing para- 
metric correlations is rather natural within the formal 
framework presented in the subsequent Section. How- 
ever, it is less obvious on an intuitive level, especially so 
since the reactance matrix is used much less frequently 
in the scattering theory proper than its 'cousin' the T- 
matrix (a very detailed exposition of the subject can be 
found, however, in 0]). Nevertheless, it can be quali- 
tatively understood based precisely on the fact that re- 
actance matrix is the Hermitian analog of the T-matrix. 
Indeed, the latter is conventionally defined by the equa- 
tion V^lxout) = ^Ixo), where |xo) is the incident plane 
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wave, and |xout) is the scattered wave, normalized to unit 
flux. Similarly, TZ satisfies an equation of identical struc- 
ture, T^|<y5v) = TZ\(Pq), where \(pv) and |(/?o) are now the 
conventionally normalized respectively, perturbed, and 
unperturbed, states in a closed system. Since we con- 
sider the parametric correlations of the real eigenvalues 
of Hermitian Hamiltonians in closed systems, reactance 
matrix is indeed quite a natural object to appear. 

The eigenvalues of TZ are expressed as tangents of the 
corresponding eigenphases (scattering phase shifts) of S. 
The parameters ya introduced in Eq. H12|l thus have an in- 
tuitive physical interpretation as the corresponding phase 
shifts. The expression of Xa (or ya) in terms of the re- 
actance matrix should be regarded as the analog of the 
microscopic definition of xq given by Eq. |(SJ). Moreover, 
the fact that the 'natural' parameters are the eigenval- 
ues of the corresponding reactance matrix allows for a 
more direct definition of the special set of curves in the 
space of Hamiltonians which were discussed above. These 
curves can now be represented as sequences of Hamilto- 
nians which map onto commuting sequences of reactance 
matrices. 

It is possible to complete the analogy with the case 
of extended perturbations by establishing a counterpart 
of Eq. 0). The distribution P{u) of level velocities u — 
dei/d\y\ can be extracted from Eq. with the help of 
the identity 

P{u) = lim |y|i?ii(f^/A = u\y\,y) . (14) 

Making use of this identity, one obtains the following 
result 

poo 

P{u) = Re / dC. 6*^"+'^ e^S»'"(i+*^'^") , (15) 

where Ka = ya/|y|- Thus, the cumulants Cm of the level 
velocity are given by 

By choosing |y| as the variable, we have already effec- 
tively performed a rescaling analogous to (up to a 
factor of tt), with C2 = C{0). Any set of additional r — 1 
cumulants is sufficient to invert Eq. (|16|l . thus provid- 
ing the remaining (non-linear) variable transformations 
between phenomenological constants Cm and universal 
parameters ya. Equation H16() establishes the universal- 
ity of Eq. ©. 

Although we have kept the discussion above rather 
general, the most interesting application from the prac- 
tical point view is when V is diagonal in the coordinate 
representation, and Hq and H^is together describe weakly 
disordered metals. The latter are characterized by the 
hierarchy of scales L £ ^ Xp. In this case it is conve- 
nient to distinguish two subtypes of local perturbations 
V. In the first subtype, the spatial extent ly of V is 



smaller than £. For such perturbing potentials, Eq. ® 
can be more compactly rewritten as 

crioc(A, Ai; x) = tr (1 + iUXi) (1 -I- iTZXy^ , (17) 

where TZ is the reactance matrix projected onto the Fermi 
surface, and tr denotes the trace in the space of the scat- 
tering channels. The full reactance matrix is defined by 
the equation 

fR== (l-FRe(G"'))-V , 

where the retarded average Green function is defined as 
(G"*) = ((e_F - Ha - i?dis + «0)-i), so that d\ is for- 
mally distinct from the average reactance matrix ((1 — 
V'R,eC^)~^V); nevertheless the difference between d\ and 
the average reactance matrix enters only in the sublead- 
ing order in Xp/t, and can be ignored. Transforming to 
the momentum representation it is convenient to split off 
the angular variables in £H as 5H(p, p') = 5R(^, n; n'), 
where f — p^/2m ^ ep- The projected reactance matrix 
TZ is defined as TZf^^ii' = 7rf'fH(0, n; 0, n'). 

In the opposite hmit L ^ ly ^ £, Eq. (fTTjl retains its 
structure provided TZ is replaced with 

TZp = P'/\Om,i^;e,n')P^^\e), 

where 

e + (i/2t)2 ' 

T = £/vp is the mean free time, and the trace opera- 
tion tr is redefined to include the ^ degrees of freedom. 
In other words, the projection onto the Fermi surface is 
smeared over a strip of width 1/2t around it, reflecting 
the fact that the average Green function exponentially 
decays over distances larger than £ so that a scatterer 
of size larger than £ is effectively split into independent 
chunks of size £. Since the algebraic structure of the 
parametric correlation functions is identical in both cases 
ly < £ and ly > £, below we will often omit the subscript 
P. 

In summary, the two ways of defining y either in terms 
of the solutions of Eq. (|16|l or as the eigenphases of the 
reactance matrix stand in direct correspondence with the 
definition of xo either phenomenologically Q or micro- 
scopically |(SJ|. While very convenient for theoretical anal- 
ysis, the definition in terms of the eigenphases is of little 
practical utility since measuring TZ is hardly feasible. On 
the other hand, cumulants of the level velocity distri- 
bution are readily accessible in appropriately designed 
experiments. 

The universality of the parametric correlation func- 
tions is explicitly demonstrated through the existence of 
their translationally invariant parametrization in terms 
of the scattering phase shifts. In this context, the stan- 
dard universal Wigner-Dyson statistics of energy levels 
{spectral points) could be viewed as the boundary condi- 
tion to the universal statistics of spectral lines as func- 
tions of coordinates in the space of Hamiltonians. 
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A. Orthogonal ensemble 



and H17|l are replaced with 



All of the preceding discussion can be applied verbatim 
to the case of orthogonal symmetry provided Eqs. Q) 



J 



= l + Rey d\id\2 j dA J(A,Ai,A2)e*'^'*+(^i-^^-^)-^*°'(^'^i'^^'^°'^) , (18) 

I 



where 



J(A,Ai,A2) 



(AiA2-A)2(l-A^) 
(A2 + A2 + A2 - 2AA1A2 - If 



and the global and local parts of cr^"-' are given by, re- 
spectively, 



(A, Ai, A2; X,) = ^ (1 + 2A?A2 - A^ - A? - A^) 



and 



af:)(A,Ai,A2;x) 

= itrlni±^^^^i^^^^S±^. (19) 

As before, 7?. has to be replaced with TZp with the corre- 
sponding redefinition of tr if ly S> I- 

A similar set of expressions can also be written in the 
symplectic (/3 — 4) case. 



B. Level velocity distribution and Berry's 
conjecture 

Applying Eq. lO to Eq. lO, we find that Eq. is 
a special case of a more general expression: 



P(u) = ^Re / dC e 
2 Jo 



J7r/3u+C/2 g-^(/3/2)i:Mn(l+jCK„) 



(20) 

It is interesting to note that the distribution of level ve- 
locities can be obtained using a completely different ap- 
proach. As can be easily seen from the expansion of an 
arbitrary eigenvalue of iJ + XV to linear order in X, 

V,, ^ [ dr |V'.(r)|2T/(r) , 



the level velocity coincides with the expectation value of 
V in the i-th eigenstate of H. Restricting our attention 
for simplicity to the case ly ^ £, we can employ Berry's 
conjecture about the distribution of wave function values 
in chaotic systems [T9ll20) : 



1 



(det/-i)/5/2 



,-(/3/2)Jdr<ir>t(r)[/-i]„,V(r') 



where in d spatial dimensions the matrix elements of the 
integral operator / are given by the Friedel function 111(1 , 
which in the limit <C ^ can be approximated as /d(r — 
r') = J dn e*P^^" (''""" ) — Jd/2-iiPFr)- In this expression 
the integration is performed over all directions of the unit 
vector n and r = |r — r'|. The fields are complex in 
the unitary case, and real for f3 = 1. While originally 
conjectured for chaotic systems '19], this expression was 
recently proved 20] to hold locally in diffusive systems 
of unitary symmetry. We thus immediately find 

P{u) = AY,{S{u - Vu/^)5{e - e,)> 



27r 

°° dC 



00 



27r [det(l + 2</l///3A)]/3/2 



(21) 



The last expression can be shown to coincide with 
Eq. (|20|) , thus indirectly confirming the validity of Berry's 
conjecture in diffusive systems of orthogonal symmetry. 
Indeed, ya/\y\ coincide with the eigenvalues of the Fermi 
surface projection oi dlZ/ dX\x^o = V. The correspond- 
ing projected matrix elements are / dr e*P^("~" '>''^V{r). 
Using the integral decomposition of /d(r — r') and the 
cyclic invariance of the determinant, the equivalence of 
Eqs. H20|) and H21|l is immediately established. This anal- 
ysis can also be in an obvious way extended to the case 
L > /y > 



C. Examples 

Before turning to derivation of the above results, it is 
instructive to analyze several special cases. For example, 
if V models a collection of several impurities, it can be 
approximated as 



nr)--E 

k=l 



VkS{r - Tk) 



(22) 



Crudely speaking, a bistable impurity would have r = 
2 and vi — ~V2, while a local defect created, e.g., by 
irradiation would correspond to r = 1. 



s 



X 



FIG. 3: DoS-DoS correlation function for tlie rank r — 1 prob- 
lem, Rii{s,y), for different strengths tany of the perturbing 
potential: tany = 0.2 (dashed line), tany — 0.8 (dotted line), 
tan 2/ = 4 (dot-dashed line), tany — > oo (solid line), while 
the bold line corresponds to the GUE Wigner-Dyson result, 
y = 0. The corresponding Fourier transforms are shown in 
the inset. 



Throughout this subsection we largely restrict the dis- 
cussion to the unitary ensemble. The corresponding re- 
sults for (3=1 can be straightforwardly obtained in terms 
of the corresponding A integrals starting from Eq. H19|l . 

Despite the fact that TZ (or TZp) is a continuous in- 
tegral operator, the structure of Eq. H22|l ensures that 
it possesses exactly r non-zero eigenvalues. Using the 
cyclic invariance of the trace, the latter are easily seen 
to coincide with the eigenvalues of the asymmetric ma- 
trix Mkk" = X]fc'=i^(rfc,r/c')/d(rfc' - rfc")i where, for 
arbitrary ly ^ L, the Friedel function /(i(r) is given by 
Eq. ini). 

In the simplest r = 1 case, Mkk' is a number, tany = 
v/{l — av), where for brevity we denote a = Ke{G^{r = 
0)) ///, and the correlation function in the unitary ensem- 
ble can be explicitly written in the form (assuming y > 0) 



Riiis.y) - 1 



d_ 

ds 



-ITS cot y 



eis) 



kis)] 



ds'e^'''™'^fc(s') 



(23) 



where the function k{s) has been defined in there- 
fore reproducing the result obtained previously .21] using 
RMT methods. The function Ru for r = 1 is plotted in 
Fig. O for several different values of y together with its 
Fourier transform. One can clearly see the gradual broad- 
ening of the central peak inherited from the (5-function 
term in Eq. jSJ . At small values of y the dominant con- 
tribution to the broadened peak comes from the corre- 
lations between a level and its parametric 'descendant'. 
It is also worth noting that the parametric correlation 



FIG. 4: Dependence of the energy levels on the overall 
strength X of an external local (i.e., localized on a single site) 
perturbation V for a time-reversal invariant tight-binding 
20 X 20 lattice model (arbitrary units). Level anti-crossing 
and the property et < £i{X) < ei+i l2j| are clearly shown. 



function inherits from its non-parametric limit the sharp 
oscillatory behavior which is reflected in the singularity 
(cusp) in its Fourier transform 
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/oo 
ds e-2--*i?,,(s,y) 
-OO 

1 + tan2 yg^{t) 



S{t) + mm{t, 1) + - In 



l + tan2y(l + 2t)2j ' 



where g{t) = 2t — 2min(t, 1) + 1. Although not directly 
obvious from Eq. (|23ll . the perturbed levels in the r = 1 
case possess the property (for y > 0) Si < Si{y) < Ei+i 
PH . This feature is clearly illustrated in Fig.|^ 

In the r = 2 case we restrict our attention to modeling 
a bistable impurity, thus setting vi = —V2 = v > 0, and 
ri2 ^ £. Denoting as before a = Re(G"(r = 0))/P, and 
also 7 = Re{C{r2i))/P, rkk' = \rk - r*;'!, and fkk' = 
fd{rkk'), we find 



tan?/i,2 



1 — ^2(72 — q;2) 



In the limit ri2 both eigenvalues, as expected, 
vanish. We omit somewhat lengthy explicit expressions 
for Rii and its Fourier transform, presenting instead in 
Fig. |S1 the corresponding graphs at several different val- 
ues of y. The Hmit can also be used to extract 
the distribution of level velocities due to a moving im- 
purity. Expanding 1 — /12 ~ (l/2d)p|,r^ wc find (set- 
ting for simplicity a = 7 = 0) that the distribution of 
u = dsi/ d{pFT') is Poissonian: 



Pin) = — e 



\^/d/v 



(24) 
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FIG. 5: DoS-DoS correlation function for the rank r = 2 
problem, Rii{s,y), for different values of tanj/ = tanj/i — 
— tanj/2- The bold line corresponds to the GUE Wigner- 
Dyson result, tany = 0, while tany = 0.2 (dashed line), 
tany = 0.4 (dotted line), tanj/ = 1 (dot-dashed line), 
tany — > ex (solid line). The corresponding Fourier transforms 
are shown inset. 



FIG. 6: Graph of Pstm{u/u) in the unitary ensemble obtained 
for pfp ~ 1.5 and 4.0 together with the limiting cases of Pois- 
son {pfP ^ 0) and Gaussian {ppp oo) distributions. The 
inset shows the perturbation V{r) arising from the displace- 
ment of the disk-shaped potential. 



This result can be generalized to study the response 
of the energy levels to a shift in the position of an ex- 
tended defect. An experimentally relevant application is 
the lateral motion of an STM tip over a disordered two- 
dimensional electron gas. Another example is a metallic 
scatterer inside a microwave 'billiard' such as those stud- 
ied in Ref. 

Approximating the potential U of an STM tip as a 
flat disk of radius p, U{r) = Uo6{p'^ — r^), and denoting 
the displacement of the center of the disk as A, the dif- 
ference between the potentials produced at two adjacent 
positions of the disk is given, to the first order in A, by 
V{r) = V{r, (j)) = 2UoXpS{r^ — p^) cos cf), where the direc- 
tion of the displacement corresponds to = 0. Therefore, 
neglecting the 0{Uq) corrections to /rf(r) = Jo{pFr), one 
finds that the operator fV/A appearing in the expression 
for the velocity distribution H21|) reduces to a continuous 
integral operator defined on a circle < cj) < 2tt: 

k[c^,c^') = -^Jo(2p;.psin((0-(/.')/2))cos0' . 

Since the corresponding eigenvalues Kq are symmetric 
with respect to zero, the distribution of level velocities 
can be written as 



^ .C("/fi)-(/3/2) i:,^>0 ln[l + (2«;„C//30)'] 
2^ 



where the average velocity u is given by 



/3 



1 



d0 J^{2pFP sin 0)cos{20) . 



In the unitary case, (3 — 2, the distribution function 
(Fig. EI) shows a crossover from Poisson to Gaussian be- 
havior as PfP increases, while in the orthogonal case 
(Fig. 0) the crossover is between the limiting behavior 
described by the modified Bessel function {l/7r)Ko{\u\) 
23] and the Gaussian limit. These crossovers explicitly 
illustrate the distinction between local and global per- 
turbations, and show that the global regime is achieved 
when the central limit theorem comes into force due 
to a large number of distinct eigenvalues Ka- Such a 
crossover has been observed in experiments on microwave 
resonators 2^] where the appropriate ensemble is orthog- 
onal. 

A new feature appearing in this example is that a lo- 
cal potential is described by a formally infinite number of 
the eigenvalues Ka of an integral operator. However, the 
finite extent of V guarantees that all but a finite num- 
ber of these eigenvalues are vanishingly small, so that, 
depending on the required degree of accuracy, there can 
always be defined an appropriate finite value of r. 



III. FIELD THEORY OF PARAMETRIC 
CORRELATIONS 

We turn now to the derivation of the results presented 
in the preceding Section. To explore the influence of a lo- 
cal potential perturbation on the ensemble average prop- 
erties of the weakly disordered system, we will employ a 
conventional approach based on the supersymmetric field 
theoretic formulation. Since this approach has been re- 
viewed extensively in the literature we will keep our 
discussion here concise, paying particular attention only 
to the idiosyncrasies of the present scheme. 
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FIG. 7: Graph of PsTuiu/u) in the orthogonal ensemble ob- 
tained for pfP ~ 1.8 and ppp = 5 together with the limiting 
cases of ppp and ppp oo. The inset shows u = (u^)^''^ 
as a function of pfP- 

Focusing on the weakly disordered system, our starting 
point is the set of single-particle Hamiltonians 

where ^p_eA = (p — eA)^/2m — sp is the kinetic energy 
operator a free-particle system subject to an external vec- 
tor potential A. The impurity potential U{v) is drawn 
at random from a Gaussian (5-correlated distribution with 
zero mean, and variance 

{U{v)U{v')) = -^5{v-v'), (25) 
2,t:vt 

where t is the associated elastic mean free time, r = 
I 



Z[i] ^ / D{iIj\%1)) exp <^ i / dr I 



where, by choosing the fields ijj to consist of 2p copies 
of both boson (b) and fermion (f) elements, the normal- 
ization of the generating function, Z[0] = 1 is enforced. 
The factor 2 in 2p is explained by the need to generate 
both retarded (r) and advanced (a) Green functions, so 
that el = {e^(Tl'' + crg^iO) ® af ® cr^^. The further dou- 
bling of the number of components of the superfields is 
dictated by the need to properly take into account the 
soft modes associated with the time-reversal invariance 
of the Hamiltonian (hence the notation tr for the corre- 



ijvp. In the following, we will limit our considerations 
to the diffusive regime, where the sample size L is greatly 
in excess of the mean free path i and where the wavefunc- 
tion is extended over the volume of the system. Moreover 
the energy scales are arranged in the hierarchy 

ei^ > - > -Be > A , 
r 

with Ec = D/L"^ and D = Vpr/d. In addition to the 
random potential U{r), the diffusive system is subject 
to a further (potentially short-ranged) arbitrary external 
parametric perturbation which can take n values Vi (r) . 

Our goal initially is to construct a formalism which 
allows, at least in principle (see below), to compute multi- 
point correlation functions of the DoS, 

p 

Rie;V) = (Y[i.{ef,V,)) , (26) 

i=l 

where £ and V denote p-dimensional vectors with com- 
ponents Ei and Vi, respectively, and ^{ei; Vi) is the value 
of the density of states at energy of the Hamiltonian 
i,p + U + Vi. For generahty, let us assume that p > n, 
so that some of Vi in H26|l are equal. The p energies 
in £ (which are generically different) are correspondingly 
split into n groups of size pj, each group matching a 
given value Vj. Moreover, from now on, we assume that 
the Fermi energy e p , which is included in the free Hamil- 
tonian ^p, is subtracted from the energies e^. 

According to the standard methods 0, ^2 ^ the gen- 
erating functional required to construct a field-theoretic 
representation of R{e:Y) in the case of unbroken time- 
reversal symmetry is written as a functional integral over 
an 8p-dimensional supermultiplet of complex fields: 



* + , (27) 

I 

spending subspace): 

where C — u™ (g) Efl + icr™i?||, and the i?-matrices are 
the projectors onto the corresponding parts of the super- 
space. To stay close to the conventional notation adopted 
in the literature, in the following we will denote a^^ = A. 
The conjugate superfields are defined by ■0 = ip^ L, where 
L = A (g) Efl -I- 1 (g) E22 (further details on notation 
can be found in |^). Finally, e'^ = diag(£^, . . . , ep. 



£ ^ - 4 - Uir) - Vir) 
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V = diag(Vi, . . . , Vp), and the absence of an explicit op- 
erator in any subspace always implies the corresponding 
identity operator. 

The only exception to this structure is the case p = 2, 
where it is sufficient to represent each of the two DoS 
factors by either a retarded or an advanced Green func- 
tion, with the opposite choice for the remaining factor, 
thus reducing the number of required components of the 
superfields from 16 to 8. 

An ensemble average of the generating functional Z [j] 
over realizations of the random impurity potential U{r) 
induces a quartic interaction of the fields 5*, which can 
be decoupled by means of a Hubbard-Stratonovich trans- 
formation with the introduction of 8p x Sp-component 
supermatrix fields Q(r). In the absence of symmetry 
breaking sources, the resulting action is invariant under 
pseudounitary transformations, Q i— *■ TQT~^, where T 
satisfies T^LT = L, and thus belongs to the pseudouni- 
tary supergroup U(2p, 2p|4p). Moreover, Q satisfies the 
time-reversal symmetry constraint Q = LC^ Q^CL. 

After performing the integration over the superfields 
one obtains the ensemble- averaged generating func- 
tional 

(Z[0]) ^ j DQ e-^[21 , 

where the action S[Q\ is 

•5[2] = j dv stiQ' + ^Jdr str(r|ln^-» , 

(28) 

and 

denotes the supermatrix Green function. The description 
of the structure of supermatrices and the definition of the 
supertrace (here, str ~ tre — trp) operation can be found 
in 4]. 

The action possesses an almost degenerate saddle- 
point manifold. Varying the action with respect to Q, one 
obtains the saddle-point equation 

Q{r) = ^{r\g[Q;V]\r) . (29) 
my 

This equation can be interpreted as the self-consistent 
Born approximation for the self-energy of the superma- 
trix Green function. The ambiguity involved in choos- 
ing among the different disconnected solutions of this 
equation is resolved by taking into account the analyt- 
ical properties of the average Green function: Q = A. 
In the limit e — > and y — > 0, the saddle-point solu- 
tion expands to fill the degenerate manifold generated by 
transformations Q = TAr~\ where T e U(2p, 2p|4p). 



A. The Non-Linear a Model 

In the standard scheme |J| leading to the NLcrM struc- 
ture of the field theory, the fluctuations of Q in the di- 
rection perpendicular to the saddle-point manifold are 
massive due to the inequality ppi ^ 1- Moreover, in the 
leading order in l/pp£, these fluctuations are quadratic 
and independent of the transverse fluctuations. As a re- 
sult, integration over the massive modes does not lead to 
any modifications of the saddle point action apart from 
an overall constant multiplying Z; supersymmetry fur- 
ther ensures that the constant is equal to 1. In the case 
of global parametric perturbations this scheme is pre- 
served 10, 12] due to the fact that V{r) is locally small, 
as discussed above. However, it is not a priori obvious 
that the same is true in the case of local perturbations, 
since the latter are locally large. Therefore, care has to 
be exercised in the derivation of the parametric version of 
NLcrM to demonstrate that (i) local perturbations do not 
destroy locally the saddle-point structure (i.e. that the 
distinction between the massless transverse and massive 
longitudinal fluctuations is preserved), and (ii) that any 
possible cross-couplings between massive and soft modes 
mediated by V do not modify the saddle-point action in 
the leading order in 1/ppi. 

To undertake this program, we begin by separating the 
fluctuations around the saddle-point manifold into trans- 
verse modes Q(j), which are nearly massless and slowly 
varying on the scale of the mean free path i, and massive 
modes JQm (r) , in such a way that Q(r) = Q{r)+SQm{i^)- 
The latter include the longitudinal fluctuations 6Qi as 
well as fast transverse fluctuations SQf. The need to ac- 
count for fast transverse fluctuations arises from the fact 
that local perturbations vary rapidly on the scale of the 
wavelength Xp, so that their coupling to fast transverse 
modes cannot a priori be ignored. In principle, a naive 
inclusion of fast transverse fluctuations can lead to over- 
counting, as, e. g. a 2pp Diffuson mode is a Cooperon, 
and vice versa. Nevertheless, below we will demonstrate 
that massive modes do not in fact generate any correc- 
tions to the slow mode action in the leading order in 
l/pp£, and a detailed calculation of the subleading terms 
is not needed. 

By deflnition, the longitudinal modes are orthogo- 
nal to the saddle-point manifold, therefore satisfying 
[SQi{r),Q{r)] = 0. In contrast, the slow transverse 
modes can be parametrized as Q{r) = T(r)AT~^(r). The 
corresponding free supermatrix Green function, GolQ] — 
G[£'^ = 0, SQm = Q;V = 0], obeys the relation: 

Cyo(r,r') = Re(G«(r,r')) 

- iTTDfdir - r')Q ^ j ' 

where the retarded Green function is defined as — 
{—£,p — U + 10)"^ and where /^(r — r') denotes the Friedel 
function Hll|) . 
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Separating the V^-dependent parts of the action, we 
represent it as 



S[Q]^S[Q] + SS[Q,SQr, 



where 



S[Q] = —^ / dr strQ2 + _ / dr str(r| lng-i[Q]|r) 



dr sti{r\\n(l-g[Q]v)\r) , (30) 



g[Q] = g[Q;V — 0], and the expansion 



g[Q] = g[Q] - —gmQmGiQ] 

It 



i2ry 



:g[Q]SQrng[Q]SQ„MQ] 



determines the form of 6S[Q,SQm]- In the hmit er <^ 1 
the cross-couphng between e and V can be neglected, 
so that g[Q] in the last term in can be replaced with 

GoiQ]- 

Employing the condition epT 3> 1, a gradient expan- 
sion of the first two terms in H3l)|l leads to the conven- 
tional non-linear cr-model action 

So[Q] J stT [D{VQf + iieQ] . 

Using the identity strln(l-Re (G") V) = 0, the last term 
in H30|) can be rewritten as 



Sv[Q] = - str In ( 1 + iQ7^ 



(31) 



where the operation str is assumed to include the trace 
over the scattering channels degrees of freedom. Simi- 
larly, SS[Q,SQ,n] takes the form 



SS[Q,SQ^]^~—stTSQf + 6So[Q,SQf] 
or 

+ i strln jl + -l-go[Q]SQ^go[Q]f 



1 



t:v{2t)' 



;Qo[Q]5Qmgn[Q]5Q„MQ]T^\ , (32) 



where 



r = 7^ 1 + iQTZ 



(33) 



while SSqIQ, SQf] is generated from the high-order terms 
in the gradient expansion, and can be approximated as 
sti SQ'j with a coefficient of the order of P/t. Both 
Qm and 5o[Q] in the above expressions are taken in 
the momentum representation, with the corresponding 
adjustment in the definition of str. The matrix T can 
be viewed as a supersymmetric generalization of the T- 
matrix T = {S — l)/2i, where S is the scattering matrix. 



The crucial property of Eq. H32() is that even if the 
eigenvalues of TZ grow indefinitely, T stays finite. More- 
over, the Hermiticity of TZ and Q ensures that the de- 
nominator in H33(l does not generate any singularities. 
Expanding the logarithm up to quadratic order in SQm, 
and averaging the action defined by the first two terms in 
the r.h.s. of 132|) over 5Qm, we find that the correspond- 
ing contributions to the slow mode action are small as 
l/ppf.- Consequently, 



S[Q]^So[Q] + Sv[Q] 



(34) 



represents the total slow mode action in the leading order 
in l/pp£. 

This completes the formal construction of the non- 
linear (T-model action. In the absence of an external per- 
turbation V, the functional integral is dominated by the 
coordinate-independent zero mode Qo provided Ei <C Ec 
for all i. In this limit one recovers the familiar zero- 
dimensional non- linear cr-model action 4], 



SoiQo] 



2A 



str[eQo] 



(35) 



which is well-known to reproduce the standard Wigner- 
Dyson correlation functions Q . The corresponding zero- 
mode contribution to the action describing parametric 
correlations is given by 



SviQo] 



1 



str In 1 + iTZQo 



(36) 



As will be shown in the next subsection, Eq. (|36|l repre- 
sents the dominant contribution to the action when V is 
a local perturbation. In the opposite case, the 1/g cor- 
rections to the action can compete with the zero-mode 
contribution, and their relative importance depends on 
the spectral composition of V. 



B. Local versus Global Perturbations 

While the zero-dimensional non-linear tr-model H35f) 
combined with the interaction action l|36|) represents the 
leading (zeroth) order term in an expansion in the inverse 
dimensionless conductance 1/g, terms of the next order 
may, under certain (and quite typical) circumstances, 
produce a contribution which can compete with, or even 
dominate, the contribution from 5*1/ [Qo]- The presence 
of such terms is best understood as resulting from the 
spatial deformation of the zero mode induced by the spa- 
tially inhomogeneous potential V{r). 

Formally such terms could have been accounted for 
by seeking spatially inhomogeneous solutions of Eq. 
at finite values of V. A simpler computational scheme, 
however, is made possible by the fact that, to the leading 
order, V couples linearly to the inhomogeneous modes of 
Q. Thus, the contribution of the inhomogeneous saddle 
point is equivalent to the result of a Gaussian integration 
around the homogeneous saddle point. 
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Employing a scheme which was originally introdu- 
ced by Kravtsov and Mirlin [2^ to explore the im- 
pact of higher mode corrections on the universal non- 
perturbative random matrix correlations in a disordered 
metallic sample, let us parametrize the variation of Q(j) 
on the non-linear manifold by setting 

g(r) = Toe^W/^y^e-^W/^To-i , 

where the W{r) is constrained by the condition that its 
anticommutator with A vanishes, {W^(r),A} = 0. To 
avoid counting the zero mode contribution twice, the gen- 
erators of the non-uniform transverse fluctuations W{r) 
are subject to the additional constraint / dr W{r) — 0. 
The spatially uniform rotations Tq parametrize the zero 
mode as Qq ~ TqATq"^. The effective zero-mode action 
is obtained by integrating over W. 

An expansion of the action in the powers of W gen- 
erates the hierarchy of non-universal corrections to the 
zero-mode action organized as a power series in the in- 
verse dimensionless conductance 1/g. Since our interest 
is in establishing the leading contribution to the zero- 
mode action rather than the investigation of 1/g correc- 
tions to the leading result, it is sufhcient to keep only the 
linear terms in the expansion of the action in W. 

Expanding the action up to the linear order in W, we 
obtain S[Q] ~ So[Qo] + Sv[Qo] + S'[W, Qo], with 



S'[W,Qo] 



— dr str 
2 



J dr str [-D{\/W f] 

AT,-\r\f\r)nW{r) 



(37) 



Since we have chosen to employ the real space representa- 
tion of the generators W of the non-uniform fluctuations, 
Eq. H37|) involves real space matrix elements of the opera- 
tor TZ, necessitating a switch from the momentum (scat- 
tering channels) representation employed in Eq. H31I) . 

In the case of local perturbations, integrating over W 
leads to a contribution which is simply a l/i? correction 
to So[Qo], and thus can be ignored in the present study. 
It is worth noting, however, that the presence of 1/g cor- 
rection to the parametric correlation functions stands in 
marked contrast to the non-parametric case where the 
leading corrections start at the order . 

Concentrating for the moment on the global perturba- 
tions, we note that, as discussed above, local values of a 
global V are necessarily small. Therefore, it is sufficient 
to approximate T ~ pi/2]/pi/2 -phe non-universal con- 
tribution to the zero- mode action Sw[Qo] is defined as 

'5'tv[Qo] = -ln(e"^ )w ■ 
Utilizing the contraction rule [23. [25l | 

(str[A(ri)W^(ri)]str[B(r2)M^(r2)]) 

= 2n(ri - ra) str[A(ri)P(r2) - AA(ri)AB(r2) 
- A{ri)LC'^B'^{r2)CL + A{ri)ALC'^ B'^ {r2)CLA] , 



where A and B are arbitrary supermatrices, and making 
use of the identity Tq — LC~^Tq^ CL, one obtains 



'5'w[Qo] 



(7rP)2 



1 {ttD)^ 



drdr' n(r-r') str V{r)QoV{r')Qo 



9 2tt 



E(^«tr(FqgoVLqQo 



q#0 



where the diffusion propagator is defined as 

1 v-^ e'l '' 
n I- - — y T-FT2 ■ 



(38) 



(39) 



Here, as above, g — Ec/A represents the dimension- 
less conductance of the disordered system. The sums 
in Eqs. (|2H1) and (jSHI) are restricted to |q| < l/£. 
The sum 



S[Qo] = SolQo] + Sv[Qo] + Sw[Qo] 



(40) 



represents the total zero-mode action describing para- 
metric correlations induced by an arbitrary external per- 
turbation. S'h'[Qo] involves only the global part of V, 
while Sv [Qo] contains contributions from both global and 
local parts of a generic V. Any corrections to this expres- 
sion involve terms of higher orders in g~^. 

We are now in position to analyze the relative impor- 
tance of various terms in S'[(3o]- Concentrating first on 
the case of global perturbations, we have to compare 
Eq. (|38|l to the term coming from the expansion of the 
action S'y[Qo] to quadratic order in the ^(r). (We as- 
sume that J dr Vi{r) = for all i, since a nonzero value 
of the integral can be accommodated by a simple shift 
of the corresponding frequency e^.) The corresponding 
contribution is (setting Re(G'") = for simplicity) 



S^v^iQo] = ^ / drdr'f!{r-r')stv V{r)QoV{r')Qo 



(41) 

If V{r) varies slowly on the length scale of the mean free 
path, making use of the identity / dr /|(r) = t/ttD, one 
obtains 



[go] = ^A^^ J2 {%QoV-ciQo 



(42) 



q#0 



Then, using 1/g = Ard(L/£)^ > Ar, we see that in 

(2) 

this limit the contribution from Sw dominates over Sy 
since each term in the sum over q in Sw is larger than 
the corresponding term in Sy'' by a factor 1/q^^^. The 
latter is large over the whole range of the summation 
over q. Taken together with the contribution 5*0 [Qo], in 
the zero-dimensional limit s < Ec these results recover 
the standard universal parametric correlation functions 
reported in the literature [lol ll^ . Specifically, in the 
simplest case of the two-point correlation function 
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evaluation of Eq. H42() leads to the following value of the 
coefficient C(0) in igj: 



C(0) = 8d' 



' J dvdv' V{Y)W{r 



y')V{y) . (43) 



Now, by contrast, let us consider a potential V^(r) 
which has a structure at scales shorter than the mean 
free path. The summation in l|38|l still extends only up 

to g ~ l/£, while Sy includes equally all harmonics of 
V . Depending on the spectral composition of V , either 
of these terms may be the dominant one. Crucially, how- 

(2) 

ever, both Sy/ and Sy have identical functional forms. 
Again using the case of the two-point correlation function 
as a representative example, we find that the contribu- 
tion from fast modes of V does not affect either the form 
of Eq. ||SJ| or the validity of rescaling Q . Its only effect 
is to redefine the coefficient C(0), or, equivalently, 37q , clS 
described by Eq. ((Tn|l . 



C. Local Perturbations 

The structure of the formalism developed above is 
quite general, allowing definitive statements about uni- 
versality of arbitrary many-point correlation functions to 
be inferred. However, its utility as a generating func- 
tional is at present largely limited to two-point averages: 
The parametrization of supermatrices Qq of high rank 
presents significant technical difficulties. Deferring the 
discussion of universality of many-point correlation func- 
tions to the concluding section of the paper, let us explore 
the particular case of the two-point parametric correla- 
tion function of the DoS ®. We set V = diag(yi, V2) 
instead of the more restrictive notation V = diag(0, V) 
employed in Eq. © , in order to better illustrate the uni- 
versality of the results. 

As discussed above, in the case of local perturba- 
tions the dominant contribution to the action comes 
from S'y[(5o]- Remarkably, the supertrace in 6*1/ [Qo] can 
be evaluated explicitly employing Efetov's parametriza- 
tion H: 



u 
V 



cos t/ I sm ( 





v-^ 



j44) 

where u and v are unitary 4x4 supermatrices and 9 are 
matrices of commuting variables 



9I1 

^BF 



i62 



h22 











(45) 



with Q < < TT and < Oi^2 < 00. The explicit 
parametrization of u and v is not needed here because 
they commute with V. For systems belonging to the 
unitary symmetry class (i.e. where the time-reversal in- 
variance is lifted), 62 = 0- 



Using the parametrization and S'y[Qo] can 

be compactly rewritten as 



Sy [Qo] = 2 ^^"^ + iTZa^'-e'''''^ j , (46) 

where for convenience we have dropped the subscript P 
from the diagonal matrix TZ = diag(7?,i, 7?,2). It is also 
convenient to rewrite this expression in terms of the scat- 
tering matrices 5i,2 which are related to the reactance 
matrices via 

7^l,2 = «(5i,2-i)(5i,2 + a)"' . 

In terms of 

5 = diag(5i,52^) , (47) 
Eq. takes the form 



SvlQo] = -strln 



1 



where we have used strln(5 + 1) =0. Note that the 
appearance of in gZI) is a consequence of the factor 
af'' in (gni) and of the identity 5(-7^) = S^Tl). Eq. ^ 
can be rewritten as 



- str In 1 
2 



- str In I 1 - iSa^^ tan - 



(48) 

Since S is diagonal in the RA space, and af^ is off- 
diagonal, only even order terms in the series expansion 
of the second logarithm in Eq. (|48|l gives a non- vanishing 
contribution to the supertrace. As a result, this term can 
be rewritten as 



- str In ( 1 + S12 tan'^ - 
4 \ 2 



(49) 



where S12 = diag(5nS2, iSjiSi). A similar transformation 
brings the first logarithm in Eq. (|48() to the form 



- str In ( 1 + cos ( 



Combining equations (|49|l and (|50|) we find 



SviQo] 



- str In a 

4 



iTZi2 cos 6 



(50) 



(51) 



where 



7^l2 = i (slSi - a) (^slSi + a 



and, utilizing the cyclic invariance of the trace, TZi2 = 
diag(7^l2,7^l2). 

Performing the trace over the bf, ra and TR indices 
in H51(l . and again using the cyclic invariance of the trace, 
one obtains the effective action: 



SviQo] = ^trln 



1 + 2i7^l2AlA2 - 7^?2(A? + - 1) 



(52) 
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where Ai_2 = cosh0i.2 and A = cos 9. In the unitary 
ensemble, the corresponding action can be inferred from 
Eq. (|^ simply by setting A2 = 1. 5'y[Qo] is thus easily 
seen to coincide with cioc , provided X(i SbTQ identified with 
the eigenvalues of 7?.i2. 

Finally, differentiating (|27() with respect to j 

and j\ generates the source term for the two-point cor- 
relation function of the form {str Qoaf^a^^)^ . Integrat- 
ing over the degrees of freedom contained in the matri- 
ces u and V, one obtains the general two-point correlator 
of DoS for the orthogonal and unitary symmetry class, 
Eqs. (CHI) and 0. 



In the absence of an explicit parametrization of the 
8px8p Q- matrices, it is not possible to perform a calcula- 
tion analogous to the one in Section llll CI to demonstrate 
that this is indeed the case. However, such a calculation 
is, in fact, not necessary. The translational invariance 
can be inferred instead from the analysis of the massive 
modes in Section Fill Al Indeed, it was shown there that 
locally strong perturbations do not lead to significant 
coupling between massive and soft modes. Consequently, 
without loss of generality one can absorb Vi into Hq and 
redefine the remaining potentials as Vi Vi — Vi. The 
corresponding reactance matrices would be automatically 
redefined as TZi ^ 0,TZi ^ TZu- 



IV. UNIVERSALITY AND CONNECTION TO 
RMT 



B. Connection to RMT analysis 



A. Universality 

Although a specific model of Hq and i/dis was used in 
the calculation, the results are valid for any disordered 
system whose spectral statistics exhibit the Wigner- 
Dyson phenomenology, since the only crucial 'ingredi- 
ents' in the derivation were the existence of a unique 
zero mode in the NLctM description, and a clear distinc- 
tion between massive and soft modes. The results are 
also valid for generic chaotic systems, although a special 
consideration may be needed to properly take into ac- 
count the contribution from the Lyapunov region to the 
non-universal terms [26l [2^ . 

The calculation in Section llll Cl underscores the univer- 
sality of the two-point parametric correlation functions in 
the case of local perturbations by demonstrating that the 
phenomenological parameters ya are the eigenphases of 
the scattering matrix iSi2 = S^Si. The importance of 
this result lies in the fact that S12 describes the scat- 
tering off the potential V2 — Vi when the unperturbed 
Hamiltonian is Hq -I- Vi -I- i/dis- Thus, all dependence on 
the 'reference' Hamiltonian Hq is excluded from the re- 
sult: parametric correlations between any Hi and H2 are 
parametrized by the eigenphases of the scattering matrix 
off the potential H2 — Hi with Hi playing the role of the 
unperturbed Hamiltonian (or, equivalently, scattering off 
Hi — H2 with H2 as the unperturbed background) . 

In the two-point case, as we have just shown, the cor- 
relation functions between two Hamiltonians H + Vi and 
H + V2 depend on a single reactance matrix TZi2. At 
the same time, the action HH1|I suitable for the calcula- 
tion of a generic p-point correlation function involving n 
different values of the perturbing potential Vi apparently 
depends on the full complement of n reactance matrices 
TZi. On the other hand, a straightforward generalization 
of the notion of translational invariance in the space of 
Hamiltonians from two-point functions discussed in the 
Introduction to the general case immediately leads to the 
conclusion that such correlation functions should depend 
only on the n — 1 'mutual' reactance matrices which can 
be chosen as, e.g., TZu, i = 2, . . . , n. 



The n = 2 case for arbitrary p (and 3 = 2) has been 
studied recently using RMT techniques [isL Il5| . In order 
to extend the RMT results to generic disordered/chaotic 
systems a phenomenological ansatz based on Berry's con- 
jecture 19:] was employed in [131. The present analysis 
affords an opportunity to achieve such an extension in a 
more rigorous way. 

We begin by noting that in the p — 2 case the struc- 
ture of the 2-point parametric correlation functions ob- 
tained above for the disordered systems with broken 
time-reversal invariance is identical to the structure of 
the results obtained in [isl [l5l |. Indeed, according to 
Refs. [is, 15] the two-point parametric correlation func- 
tion in random matrix ensembles of unitary symmetry 



i?ii(fi) = i-|p-i[fc(s)-5(s)]} mi 



Here, as before, s — ^/A, the differential operator V is 
defined as 



I? = det I 1 - D^RMT 



d 

ds J ' 

the function k{s) has been defined in Q and IHrmt is the 
random matrix version of the reactance matrix. 



V 



l-Re{G^)V 



In the standard random matrix ensembles the real part 
of the average Green function is a diagonal matrix, thus 
simplifying the structure of the reactance matrix. Using 
the Fourier transforms 



and 



k{s) - S{s) = - 



1 dX 



dAi 
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we immediately recover Eqs. ((TJ and where Xa are sis based solely on the method of orthogonal polynomials 

the eigenvalues of 91. It is intriguing that the supersym- 0, ^| . 

metric structure involving the pair of compact (A) and It was shown in [13,13 that in unitary random matrix 

non-compact ( Ai ) variables is reproduced under the guise ensembles the p-point parametric correlation function has 

of the dual pair (2), 2?^^) as a result of the RMT analy- the form 



pi 



^PiP. ({eJLi) = ^"iU tr S{e. - H) J] ^i^f -H-V)) 

i=l j'=pi+l 



det 



, (53) 



where Si — Ci/A, the indices i and j run over the range 
1, . . . and i' and j' are in the range pi + 1, . . . ,p. 
All energies Si inside the subsets [l,pi] and [pi + l,p] 
are assumed to be different, which corresponds to ne- 
glecting the ^-function terms describing self-correlations 
of levels. The generalization of Eq. (|53|l to the case when 
some energies coincide can be found in [TsL ITsj . We now 
note that, in complement to the orthogonal polynomial 
method used in Refs. 0, the parametric correla- 
tion functions in random matrix ensembles can be al- 
ternatively studied using the non-linear sigma model ap- 
proach. The resulting sigma model action has the form of 
Eq. H36|) parametrized by 91rmt- Since they arc described 
by the same action, the correlation functions have the 
same functional form irrespective of whether the averag- 
ing is performed over H drawn from an invariant distri- 
bution or over H — Hq + H^is = + U{r) where the 



distribution of U{r) is described by Eq. (|25|l . It follows 
that the p-point correlation functions in generic disor- 
dered/chaotic systems of unitary symmetry are given by 
the universal Eq. H53|l with the operator T> parametrized 
by the corresponding reactance matrixTZ. This conclusion 
about universality and parametrization extends also to 
the level-number-dependent correlation functions stud- 
ied in 0,01 since the latter are based on Eq. 
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